Spintronics in M0S2 monolayer quantum wires 
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We study analytically and numerically spin effects in M0S2 monolayer armchair quantum wires 
and quantum dots. The interplay between intrinsic and Rashba spin orbit interactions induced by 
an electric field leads to helical modes, giving rise to spin filtering in time-reversal invariant systems. 
The Rashba spin orbit interaction can also be generated by spatially varying magnetic fields. In 
this case, the system can be in a helical regime with nearly perfect spin polarization. If such a 
quantum wire is brought into proximity to an s-wave superconductor, the system can be tuned into 
a topological phase, resulting in midgap Majorana fermions localized at the wire ends. 

PACS numbers: 71.70.Ej, 85.75.-d, 73.63.Kv, 78.67.-n 



Introduction. Atomic monolayers such as graphene 
sheetall have attracted much attention over the 
years. However, the small spin orbit interac- 

tion (SOI) in graphene makes spin effects negligibly 
smalhl^HS In contrast, transition-metal dichalcogenide 
semiconductors^j^ in particular M0S2, possess giant 
values of SOI.I^ Combined with a direct band gap 
this SOI makes these materials attractive for optical 
effect s.'2SH2Sl Howe ver, p revious work emphasized val- 
leytronics in MoS2,'^^'^ while the spin degrees of free- 
dom have received much less attention, despite the fact 
that these materials can be expected to display interest- 
ing spintronics effects, such as helical states in quantum 
wires, Majorana fermions, spin qubits in quantum dots, 
electrical control of spin, etc. This gap of understanding 
has motivated the present work where we will propose 
and analyze spin effects specifically for quantum confined 
structures in M0S2 monolayers. 

One of our main findings is that the intrinsic SOI 
needs to be complemented by Rashba-like SOI to ob- 
tain interesting spin effects in the conduction band. In 
particular, we will focus on suitably defined quantum 
wires of armchair type and show that they allow for 
helical modes, with and without time-reversal symme- 
try. The Rashba-like interaction can be generated by 
breaking structure inversion symmetry with gates or 
adatoms or, alternatively, by nanomagnets with alter- 
nating magnetization direction. Helical modes serve as 
basis for spin filters^Sl but also as platform for exotic 
quantum states such as Majorana fermiona^^or fraction- 
ally charged fermions .l^Sl We finally discuss quantum dots 
with well-defined Kramers doublets that can serve as spin 
qubits.!^ 

Bandstructure. A molybdenum disulphide (M0S2) 
monolayer consists of two layers of S atoms stacked over 
each other forming an effective trigonal lattice and of 
Mo atoms located in the center of the sulphur lattice, 
see Fig. [l^. The Brillouin zone consists of a hexagon 
with two nonequivalent corners (valleys) at K (t^ — 1) 
and — K (r^ — —1) (see Fig. ITb) that determine the low 
energy spectrum of the monolayer.^'' 30 33 'pj^jg p^^t of 
the spectrum is dominated by three d orbitals of Mo: the 
conduction band by Jt/'c) = \dz^) and the valence band by 



Ionian is given by 



\dxy)) / \/2. The effective Hamil- 



Hq = hvpikxTzCTl -f kya2) + A(T3, 



(1) 



where the Pauli matrices Ui act on the d orbital space. 
The momenta k^ and ky are calculated from the cor- 
responding valley characterized by t^. Here, vp ~ 
0.53 X 10^ m/s is the Fermi velocity and the mass term, 
A w 830 meV, arising from broken inversion symmetry, 
has been extracted from DFT calcnlations.'22 The spec- 
trum of Hq is given by Ec,v = '^\ {f^'VpYikl + k^) + A^. 

The large gap 2A between the valence and conduc- 
tion bands makes the monolayer attractive for optical 
effects. 123125] ^pj^g wavefunctions at fixed energy E and 
momentum ky are written in the basis {ipc, ipl") as 



'<Pr,,pix) 



„i{T^K+pk!c)x+ikyy 



(2) 



where 6r^^p — hvpiTzpkx — iky)/{E — A), kx > 0, and 

the index p labels right- and left- movers in x direction. 

We focus now on the quasi-one-dimensional limit 

where the system forms a quantum wire. Similarly to 




FIG. 1. (a) The M0S2 monolayer lattice consists of Mo 
(large green dots) each connected to six S (small blue dots). 
The armchair quantum wire in the monolayer can be formed 
by metallic gates (yellow area) that fix the propagation direc- 
tion (defined as the y axis) to be perpendicular to one of the 
lattice translation vectors, say ai (red arrow), (b) Brillouin 
zone where the valleys K and —K lie on the fcj, axis which 
is perpendicular to the direction of propagation given by the 
ky axis. Note that the boundaries of a M0S2 flake are not 
important for this setup. 
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FIG. 2. (a) The energy spectrum of Ho [e(fca) = E{\^kya) — A] for a quantum wire of width W = 50a as obtained by numerical 
diagonalization. Parameters are chosen as t = 1.27 eV, A = 0.83 eV, and a — 0.32 nm. All levels are degenerate only in spin. 
This spectrum is in good agreement with our analytical predictions, see Eq. ([3|. (b) The Rashba SOI term Hrx, Q-r = 10 meV, 
lifts the spin-degeneracy and results in the spin-dependent shift of the wavevector ky. (c) The remaining SOI terms: intrinsic 
SOI Hso, a = 38 meV, and Rashba SOI Hay, an = 10 meV, lead to the anticrossings in the spectrum (red dashed circles). 



carbon-based niaterialsp3E51this quantum regime can be 
achieved in two ways: either by growing a nanoribbon of 
M0S2 with particular boundaries or by electrostatically 
confining the electrons into a quantum wire by placing 
metallic gates on a M0S2 monolayer flake (with unspec- 
ified boundaries), see Fig. [l] In both cases we consider 
the armchair regime where the direction of propagation is 
perpendicular to a lattice translation vector. The details 
of the boundaries are not essential provided that they do 
not suppress the transport. In addition, we assume for 
both cases hard-wall type boundary conditions. A most 
characteristic feature of such armchair quantum wires is 
that the two valleys [±K — (±47r/3a, 0)], projected onto 
the propagation direction along the y axis, coincide at 
k = 0. The valleys easily hybridize (lifting their degen- 
eracy), for instance, by impurities, irregular bound aries 
or, in particular, by our hard- wall boundaries .S^I^ni 

We determine now the spectrum for a quantum wire 
of width W = Na, where a is the lattice constant and N 
the number of unit cells in the x direction. To find the 
quantization conditions on k^, we virtually extend our 
quantum wire by two sides, W = {N + 2)a, and impose 
the boundary conditions at these virtual sites on the to- 
tal wavefunction ipix) = J2t,,p '^r,,pipT,,p{x) = J2j V'i(a:), 
where iJjj{x) is the probability amplitude to find the elec- 
tron in one of the orthogonal states j = d^^, d^^-y^, d^y 
This gives us six conditions (three at each edge) for 
four fundamental solutions [see Eq. ^]. To solve this 
overconstrained boundary value problem we use mixed 
boundary conditions. On the orbitals ^^2 and d^^-y^ we 
impose Dirichlet boundary conditions, ipd 2,d 2_ 2(2^ — 
0, W) = 0, while on the orbital dxy we impose von Neu- 
mann boundary conditions, d^ipd^ (x = 0, W) — 0. 

These boundary conditions can be fulfilled only for 
those values of k^ = \Km\ that satisfy {K + Km)W — Trm, 
where m is an integer. If W^ = (3M -I- l)a, where M is 
a positive integer, this leads to k^ = Trm/W. We note 
that in this case all energy levels except the lowest one 
in the conduction band and the highest one in the va- 
lence band are two-fold degenerate. If ly = {3M + 2)a 
[W = 3Ma], this leads to k„ = {nm + 2tt/3)/W' 
[Km = {nm — 27r/3)/Vl^']. In this case, the energy lev- 



els are non-degenerate. The conduction band spectrum 
for small momenta is quadratic. 



E — 



^{hvpHmY + {hvpkyf 



A2 



{Hvpky)"^ 

2A„ 



■(3) 

where we define the minimum energy for the mth sub- 
band as Ajn = yJ{hvFKmY + A^, sec Fig. pk. We note 
that the slope of the spectrum branches at small mo- 
menta is decreasing with increasing m, resulting in cross- 
ings of different subbands, see Fig. [2J:. The correspond- 
ing wavefunction is given by 



V'm(a;) = ipi,p{x) - ij-i-p{x), 



(4) 



where p = sguK^- Again, we note that for W = (3M -|- 
l)a the subbands m and —m have the same energy. 

To confirm our analytical results numerically, we de- 
velop a tight-binding model for a honeycomb lattice com- 
posed of two kinds of atoms, A (representing dz2 orbitals) 
and B (representing dx2_y2 +iTzd,j.y orbitals). The effec- 
tive Hamiltonian Hq consists of an on-site energy term 
and a term describing hopping between nearest neigh- 
bours, respectively. 



Ho 



/ , ^i'^i/i'^V 



7 . ''ij'^ip.'^jfJ.' 



(5) 



where c] creates an electron with spin fi at site i. The 
on-site energy Si is equal to A (—A) on A {B) sublattice. 
The hopping matrix element tij is assumed to be uniform, 
tij = t = 2hvp/\/3a. The Hamiltonians Hq and Hq are 
equivalent for momenta close to ±K and result in the 
same low energy spectrum, see Fig. [2k. 

Intrinsic spin orbit interaction. The intrinsic spin orbit 
interaction in M0S2 monolayer is much larger than in 
other monolayers, for example, in graphene, arising from 
d orbitals of the heavier atom. Symmetry arguments 
confirmed by DFT calculations lead to the intrinsic SOI 
Hamiltonian of the fornPH 



H, 



zSz{l -c^a), 



(6) 



where Pauli matrices Si act on spin space, and a — 
38 meV is the SOI strength. 

Rashba spin orbit interaction. Breaking of structure 
inversion symmetry by an electric field E along the z 
axis perpendicular to the monolayer leads to a Rashba 
term of the form^ Hj^ = Hjix + Hf 



iRy, 



Hbx = ~aRSx(J2, Hfjy = ttRT^Syai, 



(7) 



where the Rashba SOI strength, in general, is propor- 
tional to the electric field strength. Such an electric 
field co uld be produced by gates'^ or by doping with 
adatoms.l^^Ml Por both cases, api is best determined by 
ah initio calculations or experimentally. 

An alternative way to generate Rashba SOI is to ap- 
ply a spatially varying magnetic field,l22l produced, for 
instance, by nanomagnets.SSBD For example, a magnetic 
field B„ rotating in the plane of a quantum wire produces 
the Zeeman term 



m 



Az[s^ cos(fc„y) + Sy sin(fc„y)], 



(8) 



where A^ — g^sBn/'^, and the period of the rotating 
field is A„ — 27r/fc„. The unitary spin-dependent trans- 
formation [/„ = eKp{—iknySz /2) allows us to gauge away 
the coordinate dependent term H^ in the Hamiltonian 



H„ 



f iJ^. This resuhs in H, 
Ho + Hso - a_R„s^cr2 + 



Azsx, (9) 



where ann = fivphn/'i., so the strength of the induced 
Rashba SOI depends only on the rotation period A„ but 
not on the magnetic field strength _B„. We note that 
the induced Rashba SOI described by -fffln — —ajinSz<y2 
reaches ajin ~ 11 meV for the nanomagnets placed with 
a period A„ = 100 nm. Here we note that for a quantum 
wire created by gates we can estimate that the misalign- 
ment angle should be less than a/A„ (i.e. < l°).l33ilf one 
works with a nanoribbon, then the propagation in the 
armchair or zi gzag d irection should be favoured by the 
growth process.H2H2l ^fi/Q note that also a Zeeman term 



H. 



'AzSx, which breaks the time-reversal invariance 



of the system, inevitably arises. 

To account for the Rashba SOI of Eq. (l7| in the tight- 
binding model, we allow for spin-flip hoppings,'^'^^' 



Hr. 



Siaji 



E 



c^e. 



(Bij X Gz) 



^fifi'Cjti' , 



(10) 



where s = (sx, Sy, Sz), and where the unit vectors e^ 
points along z and e^j along the bond connecting sites i 
and j. The intrinsics SOI, see Eq. (l6|, can be modeled 
by 



H. 



2ia 
3n/3 



E 



^ij'^i^i^z,iJ.iJ.'Cj,_i', 



(11) 



where the sum runs over the next-nearest neighbour sites 
belonging to the B sublattice. The spin dependent am- 



plitude Vi 



--ji 



±1 depends on whether the elec- 



tron takes a right or left turn by hopping from i to j 



m 



These two terms Hso and Hr are constructed in such a 
way that they are equivalent to Hso and Hr in the low- 
energy sector. We note that by taking only part of Hr 
and changing Sx to Sz^ we can model Hun- The Zeeman 
term Hz is given by 



H, 






(12) 



Spectrum with SOI. A part of the Rashba SOI, Hr^ 
{Hr„), can be easily included in Hq. The spin s^ (s^) is 
a good quantum number for the Hamiltonian Hq + Hrx 
(Ho + HRn). In this case, the SOI only results in the spin- 
dependent shift of the momentum, ky ^ ky — SxOlr/Hvf 



{ky 



SzCtRn/fiVF), see Fig. pb. Next, we treat the 



remaining SOI terms, H' = Hso + HRy (Hso), as a per- 
turbation. We note that H' (Hso) is proportional to r^. 
At the same time, the wavefunctions t/jmix) [see Eq. Q] 
are eigenstates of the Pauli matrix ti, so the intrasub- 
band matrix elements vanish,S31 which is consistent with 
Kramers degeneracy at ky = for a time-reversal in- 
variant Hamiltonian. The intersubband matrix elements 
imm' , however, are non-zero, but they contain a strong 
suppression factor arising from the sublattice degree of 
freedom as follows. At small momenta, the mass term 
A(73 dominates in the Hamiltonian, so the wavefunctions 
^Jm{x) are close to the eigenstate of 0-3. As a result, 
the sublattice terms in H' (Hgo), 1 — era and cti, lead to a 
suppression of SOI effects, where the intrinsic SOI is sup- 
pressed by a factor {E — A) /A <C 1 and the Rashba SOI 
by a factor \/{E — A)/A. Thus, the corrections to the 
spectrum are small in the parameter tmm' /^mm' , where 
^mm' denotes the subband splitting at given ky. How- 
ever, these terms lead to an anticrossing between two dif- 
ferent subbands with opposite spin (with the same spin) 
along X (along z), see Fig. [2J:. We note here that in spite 
of having strong SOI, the spin degeneracy is not lifted in 
case of quantum wires. 

Helical modes via electric field. The Rashba SOI in- 
duced by an electric field offers the possibility to gen- 
erate helical modes in a time-reversal invariant system. 
As shown above, H' results in subband anticrossings, see 
Fig. [2];. Sufficiently far away from them, tmm' ^ '^mm' : 
the subbands are spin-polarized by the Rashba SOI Hr^ 
in the x direction, see Fig. [3] However, passing through 
the anticrossing the spin polarization goes through zero 
and changes sign. All this suggests that if the Fermi level 
is tuned close to the anticrossing (see Fig. l2fc) in such 
a way that there are four propagating modes (two left 
and two right), the system is in a quasi-helical regime. 
The lowest subband n = 1 is almost fully spin-polarized 
and transports opposite spins into opposite directions, 
whereas the next subband n — 2 \s only partially po- 
larized. This means that scattering due to impurities 
between subbands is allowed and helical modes are not 
protected from backscattering. 

Helical modes via magnetic field. If a Rashba SOI 
(along x) is generated by a spatially varying magnetic 
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FIG. 3. The spin polarization {s^) along x direction as func- 
tion of momentum ky for the nth level defined in Fig. |2p. The 
spin projections onto the y and z directions vanish. Away 
from the anticrossings the spin is almost perfectly aligned 
along X. If the chemical potential fi, defining the Fermi mo- 
mentum kp for the nth subband, is tuned close to the an- 
ticrossing (see Fig. [2t), the total polarization (s^) of a left 
(right) propagating electron is non-zero. 



field, the time-reversal invariance of the system is bro- 
ken, giving rise to a Zeeman term Hz- The corresponding 
magnetic field, pointing along z, is perpendicular to the 
spin quantization axis determined by the Rashba SOI. 
Thus, the spin degeneracy at ky = gets lifted, and a 
gap of size 2Az is opened, see Fig. |4] The spin polariza- 
tion along z is given by 



quantum dota^SEfll jn M0S2 can be created by gates .1^^^ 
We note that the confining potential should be sharp 
enough to lift the valley degeneracy, as shown above for 
the quantum wires, and, in addition, to insure a non- 
equidistant spectrum, which is more suitable for opti- 
cal experiments. If vanishing boundary conditions are 
imposed also along y, the momentum ky is quantized, 
Kyn = ^^^/{Wy + 2a/-\/3), where Wy is width in the y 
direction, and n is a positive integer. The dot spectrum 

then becomes -E„,m — J {hvpY{K'^^ + k^^) + A^. Each 

level is spin degenerate and the intrinsic SOI can nei- 
ther lift this Kramers degeneracy, nor, due to its symme- 
try, change substantially splittings between levels [in the 
small parameter a{E — A)/ A ^ a, see above]. The spin 
degeneracy can be lifted with a magnetic field, say, along 
z. Similar to nanotubes, - EDSR can then be achieved 
by applying an oscillatory electric field E also along z, 
causing an to oscillate and thereby inducing spin rota- 
tions at a Rabi frequency ~ \oiR\. Thus, we conclude 
that quantum dots in M0S2 host well-defined Kr amers 
doublets that can serve as platform for spin qubits.'^SESl 
We acknowledge stimulating discussions with Parisa 
Fallahi, Richard Warburton, and Dominik Zumbuhl. 
This work is supported by the Swiss NSF, NCCR 
Nanoscience, and NCCR QSIT. 



Appendix A: Majorana Fermions 



w;t 



V^it + 4A 



(13) 



1, provided that A^ < l&a]^^/A. 



where lo^^ is the energy difference between spin up and 
spin down states at given momentum ky for the unper- 
turbed problem Hq + Hun + Hso- If the chemical po- 
tential II is tuned inside the gap, there is one mode prop- 
agating to the left and one to the right. Moreover, these 
two modes carry opposite spins with almost perfect po- 
larization, \{Sz) 
Thus, the system is in a helical regime. 

Majorana fermions. Helical modes as in Fig. [4] have 
attracted considerable attention in various candidate sys- 
tems not only as a platform for spin-filters^ but also as 
a platform for generating MFs-'^Zl MFs are particles that 
are their own antiparticles. When the quantum wire is 
brought into tunnel contact with an s-wave superconduc- 
tor inducing a proximity gap Asc, states with opposite 
spins and momenta are coupled giving rise to an effective 
p-wave pairing. In the topological phase, MFs emerge as 
midgap boundstates, one localized at each end of the 
quantum wire. This phase emerges if A|- > A^^ + ^^ is 
satisfied, where /x is now counted from the middle of the 
gap. Sin ce the der ivation is similar to previously stud- 
ied cases ,12SI2II11I1I1 we defer the details to App. |a] The 
M0S2 monolayer quantum wires offer the unique possi- 
bility to probe MFs not only by transport but also by 
optical spectroscopy. 

Quantum dots. In contrast to gapless graphene,HSH2l 



We give here more details of the derivation of the 
Majorana fermions (MFs) introduced in the main text. 
Thereby we closely follow the derivation given in Ref. |3S] 
which requires a few minor modifications for the present 
case. If the chemical potential fi is tuned inside the gap 
2Az opened by the magnetic field B„ at ky — 0, the 
two propagating modes are helical, see Fig. 4. The same 
helical states can be obtained by a Rashba SOI induced 
by an electric field in the presence of a uniform magnetic 
field giving rise to a Zeeman splitting 2Az. If such a 
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FIG. 4. (a) The two lowest energy levels of Ho + Hun + 
Hso + Hz, cf. Fig. [2] The Zeeman term lifts the Kramers 
degeneracy at fcy = and opens a gap 2Az. If ^ is tuned 
inside the gap, the system is in a helical regime with a left 
(right) propagating mode with spin down (up), (b) The spin 
polarization (s^) 

lowest level. The parameters are chosen as W 
10 meV, and Az = 0.05 meV. 



50a, 



CiRn 



quantum wire is brought into tunnel contact with an s- 
wave superconductor, a superconducting proximity gap 
Age is induced in the wire. Through the pairing mech- 
anism couphng Kramers partners, the hehcal states get 
paired into a p-wave-hke superconducting state. ^^ ^■^^- 
There are no propagating modes inside the gap but there 
could exist boundstates localized at the ends of the wire. 
If a certain topological criterion is satisfied, these states 
are MFs, particles that are their own antiparticles. To 
find this criterion, we describe the system by an effec- 
tive linearized model for the exterior {x = e, states with 
momenta close to the Fermi momentum, /cg = kp) and 
the interior branches (x = *j states with momenta close 
to ki = 0).'^ The electron operator is represented as 
'^{y) = Sp=±i v=e i e^''^^^'^ px^ where the sum runs over 
the right (i?, p = 1) and left (L, p = — 1) movers and 
"^ px is an annihilation operator for the (p, x) branch of 
the spectrum. The effective Hamiltonian becomes 

H = -ihvpsxsdy + Az?/3Pi(l + X3)/2 

+ A,,772P2(1 + X3)/2 + A,,772P2(1 - X3)/2. (Al) 



in the basis 



,t 



,t 



where the Pauli matrices Xi iVi) ^-ct in the interior- 
exterior branch (electron-hole) space. 

Here, Az = gpBBn/2 is the Zeeman energy, and 
V = {dE/dhky)\k =kF is the velocity at the Fermi level. 
In the limit of strong Rashba SOI (a^n >• Az,Asc), 
the strength of the effective proximity induced supercon- 
ductivity acting on the exterior branches A^c due to the 
nearly perfect spin polarization at the Fermi wavevec- 
tor kp is equal to Asc- In the opposite limit of weak 
Rashba SOI (afl„ ^ Az), Age is getting suppressed by 
the magnetic field, A^c = Asckn/kp, where fc„ = 27r/A„ 
(kn = 2aii/hvp) for Rashba SOI induced by rotating 



magnetic fields (by electric fields). Note that the Fermi 
wavevector kp grows with magnetic field as kp ex ^ Az- 
All this together leads us to the criterion for the topo- 
logical phase given by 



Az>v^ 



p- 



(A2) 



where the chemical potential /i is now calculated from 
the middle of gap 2Az- 

Similarly, following Refs. [5S| and HSl we can also ob- 
tain the localization length of the MFs. For example, 
in the strong SOI regime and for /i = 0, the wave- 
function of the left localized MF is written in the basis 
^= (*t,*;,*t'*l) as 



^Aiiy) 



-fcL''y 



/ 



V 



I e 



ikpy '■ 
-ikpy 
g-ikjpy 
aikpy 



-k^'h 



(A3) 



/ 



^(») 



where k'^l' = {Az - A.J/fe and k^^'> = Ag^/hv^^. The 
localization length is determined by the smallest gap in 
the system ^ — max{l/fc_ , l/Zc^*^^}. Here, ^l/i-, is a cre- 
ation operator of the electron with spin up (down) , where 
the spin quantization axis is determined by the Rashba 
SOI. In the weak SOI regime deeply in the topological 
phase, the left localized MF wavefunction is written as 



*M(y) = 



/ 



V- 



te 



t/4 
iir/4 



\ 



g!7r/4 



sm{kpy)e 



-fe<='i 



(A4) 



'te 



where k^'^^ — Agc/fiv'^^. Note that both solutions ex- 
plicitly satisfy the MF condition of being self-conjugate. 
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